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Abstract 

The compressible Navier-Stokes-Poisson (NSP) system is considered in M 3 
in the present paper and the influences of the electric field of the internal 
electrostatic potential force governed by the self-consistent Poisson equation 
on the qualitative behaviors of solutions is analyzed. It is observed that the 
rotating effect of electric field affects the dispersion of fluids and reduces the 
time decay rate of solutions. Indeed, we show that the density of the NSP 
system converges to its equilibrium state at the same L 2 -rate (1 + t)~ 

4 or 

L°°-rate (1 + i)~ 3 / 2 respectively as the compressible Navier-Stokes system, 
but the momentum of the NSP system decays at the L 2 -rate (1 + t)~* or 
L°°-rate (1 + t)^ 1 respectively, which is slower than the L 2 -rate (1 + t)~* 
or L°°-rate (1 + i) _3//2 for the compressible Navier-Stokes system [17 } 121 } 16]. 
These convergence rates are also shown to be optimal for the compressible 
NSP system. 

Keywords: Compressible Navier-Stokes-Poisson, internal force, optimal decay 
rate. 



1 Introduction and main results 

The compressible Navier-Stokes-Poisson system takes the form of the Navier-Stokes 
equations coupled with the self-consistent Poisson equation and can be used to 
simulate, for instance in semiconductor devices, the transport of charged particles 
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under the electric field of electrostatic potential force [23]. In the present paper 
we consider the long time behavior of global strong solutions of the compressible 
Navier-Stokes-Poisson system in IR 3 . To begin with, we study the initial value 



problem (IVP) for the compressible Navier-Stokes-Poisson (NSP) system 

d t p + V-m = 0, (1.1) 
d tm + V- (=sam) + Vp(p) + P V$ = /iA(^) + (/i + u) V(V- (1.2) 
-A 2 A$ = p-p, lim ${x,t) -> 0, (1.3) 

p(x, 0) = p (z), m(x, 0) = m (x), x G M 3 . (1.4) 

The variables are the density p > 0, the momentum m, the velocity u — — , and 



the electrostatic potential Furthermore, p = p(p) is the pressure function. The 
viscosity coefficients satisfy p > 0, |p + v > 0, and A > is the Debye length. 
p > denotes the background doping profile, and in this paper is taken as a positive 
constant for simplicity. 

When there is no external or internal force involved, there are many results on 
the the problem of long time behavior of global smooth solutions to the compressible 
Navier-Stokes equations. For multi-dimensional Navier-Stokes equations, the H s 
global existence and time-decay rate of strong solutions are obtained in whole 
space first by Matsumura-Nishida [201 121] and the optimal LP (p > 2) decay rate 
is established by Ponce [21]. The long time decay rate of global solution in multi- 
dimensional half space or exterior domain is also investigated for the compressible 
Navier-Stokes equations by Kagei-Kobayashi [TDJ, [TT] , Kobayashi-Shibata [TB], and 
Kobayashi [T3]. Therein, the optimal L 2 time-decay rate in three dimension is 
established as 

\\(p-p,u)(t)\\ L i m <C(l + t)-i (1.5) 

with (p, 0) the constant state, under small initial perturbation in Sobolev space. 
These time-decay rates in energy space reveal the dissipative properties of the so- 
lutions, but provide no information on wave propagation. For this reason, Zeng 
[29] analyzes the Green functions of one-dimensional isentropic Navier-Stokes equa- 
tions and shows the /^-convergence to the nonlinear Burgers' diffusive wave. To 
understand the wave propagation for compressible fluids in multi-dimension, Hoff- 
Zumbrun [HIE] study the Green's function of an artificial viscosity system associated 
with the isentropic Navier-Stokes equations and derived the L°° time-decay rate of 
diffusive waves. Liu- Wang [T7] investigate carefully the properties of the Green's 
function for isentropic Navier-Stokes system and present an interesting pointwise 
convergence of solution to diffusive waves with the optimal time-decay rate in odd 
dimension where the important phenomena of the weaker Huygens' principle is 
shown due to the stronger dispersion effects in multi-dimensional odd space. This 
is generalized later to the full system later in [12] where additional new waves are 
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introduced also. To conclude, the optimal L°° time-decay rate in three dimension 
is 

\\(p-p,u)(t)\\ L ^ 3) <C(l + t)-§. (1.6) 

The long time behavior for general multi-dimensional hyperbolic-parabolic systems 
is studied in [13] for multi-dimensional case with the H s time-decay rate of solutions 
obtained, and in [TH] for one-dimensional case with L p (p > 1) time-decay rate. 

When additional (exterior or internal) potential force is taken into granted, the 
global existence of strong solution and convergence to steady state are investigated 
by Matsumura-Nishida [22]. As well-known, however, the external force does af- 
fect the long time behavior of dynamical solutions. The slower time-decay rate 
for isentropic compressible flow is investigated by Deckelnick [H [2] , which is im- 
proved later by Shibata-Tanaka [261 EZ] and Ukai- Yang- Zhao [28]. The optimal L p 
convergence rate in IR 3 is established recently by Duan-Ukai- Yang-Zhao [6] for the 
non-isentropic compressible flow as 

|| (p -p,u,9- 0oo)(*)||r*CR 3 ) ^ C ( 1 + 2 < p < 6, (1.7) 

where (p, 0, 9,^) is related to the steady-state solution, under the same smallness 
assumptions on initial perturbation and the external force. The optimal LP — L q 
convergence rate with 1 < p < 6/5 and 2 < q < 6 in R 3 was also established 
by Duan-Liu-Ukai-Yang [7] for isentropic compressible flow with external potential 
force. 

It also should be noted that for the compressible Navier-Stokes-Poisson system 
related to the dynamics of self-gravitating polytropic gas, there are also important 
progress recently on the existence of local and global weak solutions (re-normalized 
solution), the reader can refer to for instance [3], EJ [251 ED], and references therein. 

However, there are few results to our knowledge on the global classical solutions 
of Navier-Stokes-Poisson system especially for the analysis of large time behavior, 
besides the local and/or global existence of re- normalized weak solution in multi- 
dimension obtained by Donatelli [5] and Zhang- Tan [30] for Cauchy problem with 
different hypothesis on the pressure-density function. To this end, we first study 
the optimal time-decay rate of the global classical solutions to the Cauchy problem 

(H3HH}. 

To be more precisely, the main purpose in this paper is to study the existence 
and uniqueness of global classical solutions and in particular the asymptotic be- 
havior on the Cauchy problem of Navier-Stokes-Poisson system. It is observed that 
the electric field leads to the rotating phenomena in fluids motion, affects the speed 
of wave propagation and reduces the dispersion effect. This makes the Huygens' 
principle observed in [17] for compressible Navier-Stokes equations invalid here for 
the compressible compressible Navier-Stokes-Poisson system and causes the slower 
time-decay rate of the momentum (or velocity vector field) to the equilibrium state. 
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Indeed, we show in this paper that the time-decay rate of the density of the com- 
pressible Navier-Stokes-Poisson system converges to the constant state at the same 
algebraic time-decay rate in R 3 (namely, (1 +t)~i in L 2 -norm or (1 + t)~ 3 / 2 in L°°- 
norm respectively) as the compressible Navier-Stokes system, but the momentum 
decays at a slower time-rate in R 3 (namely, (1 + t)~^ in L 2 -norm or (1 + 1)^ 1 in 
L°°-norm respectively) than both the compressible Navier-Stokes system [211 [T7] 
and the compressible Navier-Stokes system with external force [2"8"j lU]. These con- 
vergence rates are also shown to be optimal for the compressible NSP system. To 
this end, we consider the linearized NSP system for density and momentum near an 
equilibrium state and investigate the spectral of the linear semigroup in terms of the 
decomposition of wave modes at lower frequency and higher frequency respectively. 
This makes it possible to analyze the influences of the rotating effect of electric field 
(caused by the internal electrostatic potential force of the self-consistent Poisson 
equation) on the qualitative behaviors of the global strong solutions, and finally to 
obtain the optimal L p (p £ [2, oo]) time decay rate of density and momentum to 
the original IVP problem (|l.l|) - (jl.4p in terms of Duhamel's principle. 

First, we have the following theorem on the global existence and large time 
behavior of classical solution to the IVP problem f lLip -f TOl) . 

Theorem 1.1 Letp'(p) > for p > 0. Assume that (p -p,m ) £ //'(R^nL^R 3 ), 
/ > 4 ; with 5 =: \\(po — p, , mo)||_H" i (R 3 )nL 1 (R 3 ) small. Then, there is a unique global 
classical solution (p, m, $) of the IVP (\l.l\i-(l~4\) satisfying 



p-pe C°(R+, H l (R 3 )) n C\R + , H l ^(R 3 )), (1.8) 
m £ C°(R + , H l {R 3 )) n C\R + , H l ' 2 (R 3 )), (1.9) 
$ £ C°(R + , L 6 (R 3 )), V$ £ C°(R+, H l+1 (R 3 )), (1.10) 



and 



L 3 fc 

\\d x (p- p)(t)\\ L 2 {R:i) < C(l + t)" 3 ~2||( Po - p,m )\\ H i iRS)nL i {R3) , (1.11) 
\\d k x m(t)\\L Hm <C(l + t)-^\\(p 

- P,"io)||f/'(iR 3 )nL 1 (iR 3 )) (1-12) 

\\d k x V<5>(t)\\ L 2 {m < C(l +t)-^-^\\(p - p,m )\\ HimnLim , (1.13) 

for k = 0,1, where C > is a positive constant independent of time. 

It should be noted that the time decay rates above are optimal. Indeed, we 
shall establish the lower bound time decay rate for the global solution. 



Theorem 1.2 Let p'(p) > for p > 0. Assume (p Q - p,m ) £ H l (R 3 ) n L^R 3 ), 
/ > 4, with 5 := || (p — P? m o) ||H ! (M 3 )nL 1 (R 3 ) small enough. Denote n =: p — p 
and assume that the Fourier transform no = ^(uq) satisfies |^o(0l > cq > for 
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< |£| <C 1 with Co > a constant. Then, the global solution (p, m, <E>) groen fry 
Theorem \l.l\ satisfies for t > t with t > a sufficiently large time that 

ci(l + t)-i < ||(p-p)(t)||L 2(K 3) <C(l + t)"t, (1.14) 

Cl (l + t)-l < ||m(t)|| i2(K3) < C(l + *)"*, (1.15) 

c x (l+t)-^ < ||V$(t)|| L2(M3) < C(l +*)"*> (1.16) 



where C\,C > are positive constants independent of time. 

With additional regularity given for the initial data, we can also prove the 
optimal LP time decay rate for the global classical solution. 

Theorem 1.3 Let p'(p) > for p > 0. Assume (p - p,m ) G if^M 3 ) H L^M 3 ), 
/ > 5, with S =: || (po — p, ^o)||_ff i (R3) n Li(K3) small enough. Then, there is a unique 
global classical solution (p,m, $) of the IVP (\1.1\) -(1.4) satisfying 



||(p-p)(*)||lp(r») < C(! +^)"^ (1_ ^ ) ||(Po -p, m o)||^(K 3 )nL 1 (K 3 )' ( L17 ) 
||^WIIlp(r3) < C(l + t)"^ (1 "p )+ 5||(p -p,m )||^ (R 3 )ni i (M3) , (1.18) 

\\Vm)\\ LP{R3) < C(l + ty^'^Upo - p,m )\\ Hl(R3)nLim , (1.19) 
for p G [2, oo] . 

Remark 1.1 It should be noted that given the same kind of initial data, the strong 
solution of the compressible Navier-Stokes equations exists globally in time and 
the momentum decays at the algebraic rate (1 + t)~i in L 2 -norm or (1 + t) _ § in 
L°°-norm respectively I2(% 0, O [77| [TP]/ , while the time- decay rate fll.lll) - fll.l3p 
in Theorem li.il and (11.171) - (j!.19p in Theorem implies that the momentum of 
the compressible Navier-Stokes-Poisson system decays at the slower time-decay rate 
(1 + t)~^ in L 2 -norm or (1 + t)^ 1 in L°°-norm. This is caused by the coupling of 
the electric field and velocity vector field through the Poisson equation, which also 
destroys the usual acoustic wave propagation for classical compressible viscous flow. 
As one can see in Theorem \1.2\ that the time decay rate (l-\-t)~* of the momentum 
in L 2 -norm is optimal. 

Remark 1.2 It is also natural to show that the L°°-time decay rate established 
in Theorem \1.3\ is optimal. To this end, it is possible to show that the difference 
(p — p, m, V$) — e tB Uo decays more faster than the corresponding rates in 
U.19\) . In general, however, it can be carried by the analysis on the fundamental 
solutions ( Green's function ) of the linearized Navier-Stokes-Poisson system in the 
framework of Liu-Wang under modification. 



Remark 1.3 We also note that, the similar properties also hold for the non- 
isentropic NSP system, which is left in the incoming paper. 
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Notations: In the following part of the paper, C > and Q > with i > 1 an 
integer denote the generic positive constant independent of time. 

The paper is arranged as follows. In section [2j we apply the spectral analysis to 
the semigroup for the linearized NSP system. We establish the L 2 time decay rate 
of the global solutions for both linearized and nonlinear NSP system in section [31 
In section HI we show the LP time decay rate. 



2 Spectral analysis of semigroup 

Let us consider the IVP problem for the linearized Navier- Stokes- Poisson system 

d t n + V-m = 0, (2.1) 
d t m + c 2 Vn + A~ 2 V(-A) _1 n - /iAm -(fM+ z/)V(V- m) = 0, (2.2) 
$ = A~ 2 (-A)-y lim $(x,f) -> 0, (2.3) 

|x|— >oo 

n(x, 0) = Uq{x) =: po(x), m(x, 0) = m (x), i6l 3 . (2.4) 

In terms of the semigroup theory for evolutional equation, the solutions (n, fh) of 
linear IVP problem (I2.ip - (I2.4I) can be expressed for U = (n,m)* as 

U t = BU, U(0) = U , t> 0, (2.5) 

which gives rise to 

U(t) = S(t)U =: e tB U , t > 0. (2.6) 

What left is to analyze the differential operator B in terms of its Fourier ex- 
pression A and show the long time properties of the semigroup S(t). Applying the 
Fourier transform to system (12. 5ft . we have 

d t U = A(£)U, U(0) = U , (2.7) 

where U(t) = U(£,t) = FU(£,t), £ = (6,6,6)* and A(£) is defined as 



-i? 
<(c 2 + A" 2 |e|- 2 ) -m 2 hxs - (/* + v)£ ® 6 



The eigenvalues of the matrix A are computed from the determinant 
det(A(t)-\L) = (A+(2/i+ I ,)|e| 2 ) 2 (A 2 +(2 / i+ I ,)|e| 2 A+( C 2 +A- 2 |er 2 )|e| 2 ) = 0, (2.9) 
which implies 

A = -^|f| 2 , (double) (2.10) 
A + = - 0* + |^)|e| 2 + iz V / 4(c 2 |e| 2 + A- 2 )-(2 yU + ^|e|4, (2.11) 
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A_ = - Oi + ^)K| 2 - |zv / 4(c 2 KI 2 + A- 2 )-(2 /U + ^|e| 4 . 
The semigroup e*" 4 is expressed as 

where the project operators P , P± can be computed as 



7 

(2.12) 
(2.13) 
(2.14) 



By a direct computation, we can verify the exact expression the Fourier transform 
G(£,t) of Green's function G(x,t) = e tB as 



--: e 



tA ___ e X+t p + e X-tp_ + e A tp o 



(2.15) 



( ~ 

V 



A + -A_ 
ig c 2( e A +«-e A -«) 



i^ t {e x + t -e x - t ) 
A + -A_ 



X+-X- 



+ 



Ait A t 
s + — e ~ 

A+-A_ 







A 2 |?| 2 



o 





g®g (A + e A + f -A_e A - t ) 



A+-A- 



(2.16) 



To derive the long-time decay rate of solutions whatever in L 2 framework or on 
point-wise estimate, we need to verify the approximation of the Green's function 
G(x,t) for both lower frequency and high frequency. In terms of the definition of 
the eigenvalues (I2.10p - fl2.12p . we are able to obtain that it holds for |£| <C 1 that 



A + e A 



A+-A_ 

A + e A + f -A_e A 
A + -A_ 



cos(6£) + (/i 



„A_i_t „A_i 



cos^-^ + H^iei 2 



-A_ 



lei « i, 



where 



lei « i, (2.17) 
lei « i, (2.i8) 

(2.19) 
|z| < 1 (2.20) 



Remark 2.4 For i/ie compressible Navier-Stokes equations, we have b = c|£| + 
(9(|f| 2 ) /or |f| -C 1. T/ras, i/ie dominant behavior of wave propagation of Navier- 
Stokes-Poisson equations at lower frequency is different from the compressible Navier- 
Stokes. The so-called Hut/gens' principle is invalid here. 

To enclose the estimates, we also need to deal with the high frequency |£| 3> 1. 
In terms of the definition of the eigenvalues fl2.10p - fl2.12p . we are able to analyze 
the eigenvalues for |£| 1. Indeed, we have the leading orders of the eigenvalues 
for |f | > 1 as 



Ao~-/i|£| 2 , 



(2.21) 
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A+ ~ - (2/i + ^)|?| 2 + + Oder 1 ), (2.22) 

A_~-^ + 0(|{|-'). (2.23) 

This approximation gives the leading order terms of the elements of Green's func- 
tion as follows 

A + eA :^- eA+t = ±e A +* [1 + e ( A - A +)*l + £e A + 4 [1 - e (A— *+)*! 

(2 24) 



■v Aii % At 
A+e + — A— e — 

A+-A_ 

_A_i_t „A_t 



~ G(l)e-^, i? > 0, |f | > 77, 
|± e A + * [! _ e (A— A + )tj _ ©(iJe-Ro*, i? Q > o, |f | > 77, (2.25) 
i e A+t [! _ e (A_- A+)t] _ 0(1)6-^*, i? > 0, |f | > 77, (2.26) 



A + -A_ 26 

where Rq and 77 are some positive constants. 

Remark 2.5 The dominant behavior of eigenvalues at higher frequency for the 
compressible Navier-Stokes-Poisson system is the same as the compressible Navier- 
Stokes. 



3 L 2 -time decay rate 

3.1 L 2 decay rate for linear semigroup 

With the help of the formula ( 12.16}) for Green's function in Fourier space and the 
asymptotical analysis on its elements, we are able to establish the L 2 time decay 
rate. Indeed, we have the L 2 -time decay rate of the global strong solution to the 
IVP problem for the linearized Navier-Stokes-Poisson system f l2.1|) - fl2.4l) as follows. 

Proposition 3.1 (£ 2 -theory) Let U = (n 0l m ) G H l (R 3 ) fl L 1 (R 3 ) ; I > A, 
and denote (n(t),rh(t)) =: U(t). Then, (n,rh,E) with E = V(— A) _1 n solves the 
IVP fl2~TD-f l2~2j) and satisfies forO<k<l that 



\\d k x (E,m)(t)\\ L2m < C(l +t)-*-H\\U \\ L i m + \\d k x U Q \ 



(3.1) 



\d"n(t)\\ L 2 m < C(l +t)~-2(\\U \\ L i m + \\d h x U \\mm), 

|l 2 (R3))- 

Proof: We are going to verify that in L 2 norm, it holds 

ll^(<)lli»CR3) = \\G(;t)U \\l HRl) = \\n(t)\\ 2 L2m + \\fh(t)\\ 2 LHR3y (3.2) 

and 



\fitU(t)\\bo® =111 ■ \ k U(;t)\\l H KD = 111 " l^-'^^H' 2 

=iii-i^(e,oiiU-) + iii-i*^(e,*)iiL> 



(3.3) 
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An easy computation together with the formula ( 12. 16ft of the Green's function 
G(£, t) gives 



n (A + e A - t -A_e A + f ) ig-m (e A + i -e A - i ) 



A + -A_ 



A + -A_ 



sin(M) r -(u+hv)\£\H 



^(l)e-^(|n | + |a |), |e|>l, 
jo(l)e-^ + ^l 2 *(|no| + |m |), |£| « 1, 
\o(l)e-^(|no| + |m |), |£| » h 



(3.4) 



with Rn > a constant here and below, and 



m(£,t) = - 



ign c 2 (e A + f -e A - t ) 
A+-A_ 



+ 



(A+^+'-A-e*-*) _ -Ao* 
A+-A_ c 



ISI 2 



+ m e 



-Ao* _ »€"0 (e A + t -e A ~ t ) 



A 2 |?| 2 A+-A- 



+ e 



-H?l 2 t 



; ^l 2t fcos(6t)-(//+|r/)^|e| 2 l 

.I„Mf|2-» 



l«| 2 

- ^ (c 2 + ^)^M e -^ + Hiei 2 * + ™ oe -^i 2 *, |e| « i, 

0(l)e~ Rot (\n \ + \fho\), |£| » 1, 
O(l)e-^l 2t (|n | + |m |) + 0(l)M e -^l 2 *, |e| « 1, 
0(l)e^(|n o | + |m o |), |£| > 1, 



(3.5) 



where and below 77 > denotes a small but fixed constant. Therefore, we have the 
L 2 -decay rate for (n, m) as 



and 



|n(*)ll 



\n(Z,t)\ 2 dt+ / |H(£,t)| 2 ^ 

<C [ e-^W\\M0\ 2 +\fh(t)\ 3 )d£ 

J\t\<ri 

+ C f e -* ot (|Ho(£)| 2 + |m (0| 2 K 

■/|CI>>7 

<C(l + t)-i||(n ,mo)||i 2(R 3 )ni] 



|m(£,t)| 2 ^ + / \rn(U)\ 2 d£ 
I€I<ij > / l€l>'7 



? / e- 2 ^(|n (O| 2 (l + |el 
'l€l<u 

-.Rot/ 



-2> 



|mo(OI 2 K 



C / e-^(|n (0| 2 + |m (0| 2 )^ 
'Kl>u 

2 



<c{i + t) 2 llK,mo)||i 2(R 3 )nL] 



(3.6) 



(3.7) 
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and the L 2 -decay rate on the derivatives of (n, m) as 

*n(t)\\h m = [ \Z\ 2k mt)\ 2 d£+ [ \Z\ 2k mt)\ 2 d£ 

J\t\<ri J\€\>ri 

<C [ e-^ + ^ 2t \e k (\n (0\ 2 + l^o(OI 2 )* 

J\t\<r, (3.8) 



+ C / e-^|e| 2fc (|n (0| 2 + |m (0| 2 )^ 

'\e\>v 



<C(l + t) § fc (||(no,m )||i 1(R 3 ) + ||(no,mo)||^ k(R 3)), 



and 



^(<)IIm 3 )= / l^(f,*)l 2 # + / \^(Z,t)\ 2 dz 

J\t\<ri J\t\>n 

<c [ e - 2 ^i 2 ^r fc (K(oi 2 (i + ir 2 ) + i^o(or)^ 

■/Kl<»? (3.9) 

+ C / e-^|er fe (|n (0| 2 + |m (0| 2 )^ 
■/|€l>»j 

<C(1 +*)-*- fe (||(no,m )||ii (R 3 ) + \\(n ,m )\\ 2 Hkm ), 

for 1 < k < i. The estimates on the is obtained via the expression ( 12. 3ft for 
< fc < / as 

ll^lliw ^ii^vc-a)-^!!!,^ = cm ■ i*- 1 ^)!^^ iQ) 

^(l+^-^dKno.mo)!!^ + ||K m ) ll 2 ^)). 
The proof of the Proposition 13.11 is completed. □ 

It should be noted that the L 2 -decay rates derived above are optimal. Indeed, 
we have 

Proposition 3.2 Let U = (n ,m ) £ H l (M. 3 ) fl L 1 (M 3 ) and assume that the 
Fourier transform = J 7 (no) satisfies |^o(OI > Q) > for < |£| <C 1 with c a 
constant. Then, as t — > +oo, the solution (n, fh,E) of the IVP (12. II) - (12.41) given 
by Proposition \3.1\ satisfies 



3 



ci(l + t)" 3 < ||n(t)|| L 2 (M3) < C(l + t)~, (3.11) 
ci(l + t)~* < \\m(t)\\ L 2 {m < C(l + (3.12) 
ci(l + t)~* < P(0IUw <C(1 + (3.13) 

where c\,C > are constants which are independent of time. 

Proof: We only deal with the estimate (I3.12p for simplicity, the argument 
applies to the others. From ( 13. 5p . we have 

fh(Z,t) = - z^o(c 2 + ^jS^e-t^^iei-* 
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sin ( bt ) r -(» 
A 2 |£| 2 b 



l?l 2 



' -i£n sm(bt) c _ (u+ y m 2 t 



(3.14) 



+ O(l)e-^ 2 \\n \ + \m \) ee 7\ + T 2 , |£| « 1, 
^ 0(l)e-^(|n o | + \fh \) = T, |e| » 1, 
with i?o > a constant and for lower frequency 

b = \VWW^>^WTWW = (a- 1 + + Oidei 4 ), lei « i. 

in terms of Taylor's expansion. 

Due to Parseval's equality ||m|| = it is enough to estimate the decay rate 
of||m(£)||. It is easy to verify 



\m(t,t)\r= / \m(U)m+ / \m(S,t)\'d£ 

'\t\<v J\Z\>v 



> 



> 



[ \rh(^t)\ 2 d^ - Ce~ 2Rot 
J\t\<v 



\t\<n 



\T 2 \ 2 d£-Ce 



-2R t 



\t\>V 



> i \\Ti\ 2 dZ-C(l + tyi -Ce~ 2Rot , 
At\<v 1 



(3.15) 



where and below n > is a small but fixed constant. By direct computation we 
have 

_ S ^(X-H + + O^H))^. (3.16) 



\Ti\'d£>C 



Applying the mean value formula we have 
sin(\-H + ^\£\ 2 t + 1 (\£\H)) 

= sin{X -H + ^\H) + [sm(X-H + ^|£| 2 i + O i mh)) - sin(\-H + ^| 2 t)] 

= Sin{x -H + ^\tft) + 2 (\Z\% 

and then 

sin 2 (A-H + ^\ 2 t + ddei 4 *)) > i sin 2 (A" 1 t + ^\ 2 t) - 3 ( [I € 1 4 t] 2 ) . 



Substituting above inequality into (I3.16P we have 
\T,\ 2 di>C 



\t\<n 



— -—^(a- 1 * + — \i\H) - o 3 mmn 

\Z\<v 1^1 2 
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\Z\<t 



— ^—^{\-H+—\i\H)d^ 



e ~{2^+u)\e,\H 



\t\<n 



-C 

=h - h- 
A direct computation gives rise to 

r -(2 M 4V)|5| 2 t 

\h\<c 

J\t\<v 1^1 
and for time t > t =: *tt with R> ^ that 



(iei 4 tm<c(i+t)-5 



(3.17) 



(3.18) 



h=c 



■ Or* 



Ac 2 



ICI<^ 2 



ie -( 2 ^)IC| 2 sin 2 (A -l t+ - | C |2 )d | C | 



R fT]t^ 



Ac 2 



sin 2 (A _i t + —r 2 )dr 



o ./ R 

R 



> 



C2(l + *)-3 / e^ (2 ^ )r2 sin 2 (A^ 1 t+— r 2 )rfr ^c 2 (l + t)"5F(t) (3.19) 
Jo 2 



with C2 > a positive constant. 

It is easy to verify that F(t) is a continuous periodic function of t with the 
period Xn. It can be also shown that there is a positive constant F min > so that 

F(t)>F min =: inf F(s) > 0. (3.20) 

Indeed, it is trivial to note that for any t > to =: there is an integer ko > so 
that t £ [A;oA7r, (A; + l)A7r]. Thus, to show (13.201) . it is sufficient to deal with F(t) 
in one time periodic domain, namely, t £ [/cgA7r, {k,Q + l)A7r] for some /co > 0. We 
deal with the case A _1 t £ [koir, k ii+ |] and A _1 t £ [/cott + f , (&o + 1) 71 "] respectively, 



and obtain for R > ^-7= that 

cv A 



inf 

te[Afc 7r+f ,A(fe +l)7r] J 



> inf 



te[Afc 7r+f ,A(fc +l)7r] \ / v^F 



V 3 77 %/77T 

+ 



Ac 2 



\/6tt 



sin 2 (A tH r )rfr 



2c-/X 2cv^ 



> < 



2c^A 



Ac 



sin 2 (A t H — —r 2 )dr, A £ £ [/cqvt, fcovr 



7i" n 



/gjr 2 2 

2cv^ 
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> 2cv ^ e *#\ sin - > 0, A i G [fc 7r + -, (fc + 1)ttJ. (3.21) 

Thus, it follows from fl3~T9l) . fl3T20|) and f[3~2Tj) that 

Ji>CiF min (l + 0~3=:c3(l + 0"3, t>t , (3.22) 

with C3 > a constant. 

Combining (13.221) . (13.151) . (I3.17P and (13.181) . we obtain the lower bound of time- 
decay rate for rh as 

H^WIIlW) = ||^(t)|U 2 (R3) > ci(l + t) - *, t > ti. (3.23) 



for some positive constants ci > and t\ > 0. 

The time-decay rate of n and E can be shown in a similar fashion. Indeed, in 
view of (13.41) we need only deal with the lower frequency for the dominating term 
g-O+i^lCI 2 * cos(6i) to obtain, after a complicated but straightforward computation, 
the lower bound of time-decay rate for n as 

\\n{t)\\ L 2 m = 11^)11^3) >ci(l + <)"*. (3.24) 



We can re-present the electric field E in terms of n and the Riesz potential, and deal 
with the dominating term — j^e - ^^'^' 2 ' cos(frt) for lower frequency to obtain the 
time-decay rate as 

\\E(t) \\ L 2 {m = \\E(t)\\ L 2 m > Cl (l + t)~i, (3.25) 
where we recall that E = — t^tl. The proof is completed. □ 

3.2 L 2 decay rate for nonlinear system 
3.2.1 Reformulation of original problem 

Let us reformulate the nonlinear system (ll.ll) - (ll.4p for (p, u) near the equilibrium 

state (p, 0) = (1,0). Denote 

n = p — p, m = m, $ = $. (3.26) 

Then, the IVP problem for (n, m) is 

d t n + V-m = 0, (3.27) 
d t m + c 2 Vn + V$ - pAm - (p + z/)V(V- m) = -/ , (3.28) 
$ = A~ 2 (-A)^ri, lim $(x,t) -> 0, (3.29) 

>oo 

n(x, 0) = no(x) =: po(x) — p, m(x, 0) = m (x), 16I 3 . (3.30) 
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where c = c(p) = yV(p) is the sound speed, and 

fo = fo(n, m, d x n, d x m, <9 2 n, dim) =: V- F(n, m, d x n, d x m), (3.31) 
F(n,m,d x n,d x m) = - A' 2 V(-A)" 1 n ® V(-A)- 1 n + ^A" 2 | V(-A)- 1 n| 2 J 3x3 
+ (p(p + n) - p(p) - c 2 n)I 3x3 

Denote 

f/=(n,m)*, f7 = Km ) i . (3.32) 
We have the equivalent form of system fl3.27p - fl3.30p in vector form 

d t U = BU + V- H, U(0) = U , (3.33) 

where the differential operator B is defined as 

B = I „ ° Q , ~ V ' | (3.34) 

y-c 2 V - cZ~ 2 V(— A) -1 -//A - (/i + z/)VV- y 1 ; 

and the nonlinear term if is expressed by 

H(U,d x U) = {p,F(U,d x U)y. (3.35) 
Thus, we can represent the solution in term of the semigroup 

U(t) = S(t)U + [ S(t- s)V - #(£/, (3.36) 

with the semigroup S(t) defined via multiplier through Fourier transformation 

S(t)U = e tB U = T- x e tA ^TXJ, A(£) = F(B)(£), £ e M 3 . (3.37) 

To establish the time decay rate of the original nonlinear problem, we need 
to decompose the Green's function G =: e tB in terms of its Fourier transform 
£?(£). Indeed, by the formula (12.161) . we can make the following decomposition for 

(n, m) = G * Uo as 

n = N ■ Uo = (A/ - + m) ■ Uo, m = M-Uo = (M + m)-Uo, (3.38) 

where 

M = ( ^1 V A; ' e-^(/ - ||) + ^^5^) 3x4 , (3-40) 
: (e A + t_ e A- t) . g 
A + -A_ V l«l 2 



= T T^f-* 0) . (3.41) 
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And we have the Fourier expression for the electric field E as 

^ = = (" |^ ® = (-^ ® AT)f> + (- ® W)8b, (3.42) 
from which we can define 

E = LU = (£ + £)U , (3.43) 

where 

C = (-^®Sf) : £=(-^®% (3.44) 

It is easy to verify that the global solution (U, E) of the IVP problem for the 
nonlinear Navier-Stokes-Poisson system fl3.27p - fl3.30p is 

U = (n,m) = S(t)U + [ S(t - r)Q{U)dr, E = V(-A)~ 1 n, 

Jo 

with Q{U) = V ■ H, which can also be decomposed as 

n = N*U + [ m{t - t) * Q(U){r)dT, (3.45) 
Jo 

m = M*U + [ M(t - t) * Q(U)(r)dr, (3.46) 
Jo 

E = L*U + [ £,{t-r)*Q{U){r)dr, (3.47) 
Jo 

By Proposition 13.11 we have the following time decay rate for linear part 

\\d%N*u {t)\\ L * < C(l +t)-f-^(||n || Ll + ||<9> |b), (3-48) 

i M 



\\d2M*uo(t)\\v < C(l + t)-*--(\\u \\ L i + WdyoWv), (3.49) 
\\d:L*u (t)\\ L 2 < C(l + t)-^(\\u \\v + ||d> |M, (3.50) 



for | a | ^ 0. Furthermore, in view of f)2.17l) - fl2.19l) and the above definition f!3.38|) 
for 9T, Xi, £, it is easy to verify for some constant c 4 > that 

m)\~0(lM\e- C4ie \ \M(0\~O(l)e~ c ^ 2t , |£(£)| ~ 0(l)e~ c ^ 2 \ |£| « 1. 

(3.51) 

Thus, applying the similar argument as in the proof of Proposition 13. 11 we are able 
to obtain after a straightforward computation (which we omit the details) that 

\\d^*u (t)\\ L 2 < C(l +t)-i(WH- J 2 1 (|| Uo || £9 + ||<9> |M, q = 1,2, (3.52) 
||^f * t*o(*)ll^ < C(i + C^'^aKlb + ||9> ||l 2 ), g = 1,2, (3.53) 
||^£*«o(t)IU=» < C(l + ty^-^-^(\\u \\ Lq + ||5>o|M, 9 = 1,2. (3.54) 
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3.2.2 Global existence and L 2 -decay rate 

We are now ready to prove Theorem 11.11 and Theorem 11.21 on the optimal time 
decay rate of solution to the IVP problem of the nonlinear Navier-Stokes-Poisson 
system f l3T271) -f l3T30l) . From (13T351) we have 

Q(U) = V ■ H(U, VU) = Qi + Q 2 + Q 3 , (3.55) 

which implies for smooth solution (p,m) satisfying ||(p — p, m)||^4 < oo that 

Q 1 :=Q 1 {U)~0{i)(d x m-m), (3.56) 
Q 2 := Q 2 (U) ~ <D(l)d x H 2 , H 2 = d x (m-n), (3.57) 
Qs := Q 3 {U) ~ 0(1)0*71 • n + 0(l)(n • £). (3.58) 

Framework of short time existence: First of all, we give the local existence the- 
ory which can be established in the framework as [2H [20], [22]. Indeed, starting 
with the equations (jl.lj) — (jl.4p . we can make use of the theorem of contracting map 
to establish the local existence. The key point is that the electric field E can be 
expressed by (11.11) . (11.31) and the Riesz potential as a nonlocal term 

E = E + V(-A _1 )g^; / mds = E + V(-A _1 )dw f puds, 

Jo Jo 

which together with the L p estimates leads to 



V(— A r )div / puds ||#fc< C || / puds 
Jo Jo 



\H k 



Then, by the standard argument of contracting map theorem as [2TI |2"U| [22], one 
can obtain the short time existence of strong solution. The details are omitted. 

To extend the short time strong solution to be a global in time solution, we 
need to establish the uniform a-priori estimates. Indeed, under some a-priori as- 
sumptions we are able to obtain the expected estimates with time decay rate for 
lower order terms, and finally enclose the a-priori assumptions. We have 

Lemma 3.3 Under the assumptions of Theorem the solution (n,m,E) with 
E = V§ of the IVP problem ( 13T27D - fl330l ) satisfies for I = 4 that 

\\d k x n{t)\\ L * < C(l + t)-M, k = 0, 1, \\d 2 x n(t)\\ L 2 < C(l + *)"*, 
\\d k x m{t)\\ L 2<C{l+t)-\~i, A; = 0,1, \\d 2 x m(t)\\ L 2 < C(l + t) _l , 

||0 x fc $(t)|| L2 < c(i + *)-i-t, k = 1,2, ||0f$(t)|| L2 < C(l + t)"t, 

and /or / = 5 that 

\\d k x n(t)\\ L2 <C(l+t)-f-f, fc = 0,l,2, ||0^(t)|| L2 <C(l+t)^, 
||0M*)IU 2 <C(l + t)-i-t, fc = 0,1,2, ||0X*)||l 2 <C(l + t)~\ 
\\%$(t)\\ L 2 <C(l + t)-i-t, A; = 1,2,3, ||0^(t) || £2 < C(l + t)" 1 , 

where C is a positive constant independent of time. 



s) 4 
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Proof: Suppose that (n,m,E) G (H 4 ) 2 x H 5 and (n,m,E) G (£T 5 ) 2 x H 6 with 
_E = V$ the electric field correspond respectively to the strong solutions of the 
compressible Navier-Stokes-Poisson system for t G [0, T] subject to initial data in 
different Sobolev space. Assume that the classical solution of the IVP (I3.27p - p.30p 
exists for t G [0,T] and denote 

Ai(*) := sup {||£>*n(s)||(l + s)i + * + ||D*m(s)||(l + s)*+i + \\D 2 x (n,m)(s)\\(l 

0<s<t, k=0,l 

+ \\E{s)\\{l + s)\ + \\{Dln,D%Dlm ) Dtm){s)\\} ) (3.59) 

and 

A 2 (t) := sup {||D*n(a)||(l + a)5 + i + \\D 3 x n(s)\\(l + a) + ||(U*n, D 5 x n, D 4 x m, D 5 x m)(s) 

0<s<t, fc=0,l,2 

+ ||££ m ( a )||(l + s )| + ||£( s )||(l + s )3 + ||D*m(«)ll(l + 

(3.60) 

We claim that it holds for any t G [0, T] that 

Ai(t) < C6, A 2 (t) < CS, (3.61) 

with 5 defined in Theorem II. H In addition, the claim (13.611) together with the 
smallness assumption on 5 are sufficient for us to prove the Lemma 13 . 31 and Theo- 
rem [TJTJ 

Next, let us prove the claim (I3.6ip . It is sufficient to show the first one in 
(13.611) since the arguments can be also applied to second one in (13.611) under minor 
modifications. The proof of the first one in (13.611) consists of following three steps. 

Step 1: The basic energy estimates. Starting with ( 13. 45ft . Proposition ^. 11 ( I3.48p . 
(13.521) . and the a-priori assumption (13.591) . we have after a complicated but straight- 
forward computation that 

\\(n-N*U ){t)\\ < f \\m(t-r)*Q(U)(r)\\dr 
Jo 

< C [\l+t- r)-*-*(||g(C0(r)|| + \\Q{U){r)\W)dr 



o 



<C / (l+t-r)-»-5(A 1 (t)) 2 (l + r)- 1 dr 
Jo 

KCil + t)- 1 ^)) 2 , (3.62) 



where we have made use of ( 13.591) and (1 3 . 5 6 1) ( T3 . 5 8 1) to estimate the right hand side 
terms as 

(\\Q(U)\\ + WQWWv) < {\\(Qi + Q2 + 0s)||} + {\\(Qi + Q2 + Q 3 )\\ia] 

< C{||Z)m||||m||ioo + ||D 2 m|| + ll m IU°°ll-D 2r? 'll 
+ ||n||||£'||ioo + ||Z)n||||n||£oo} 
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+ C{||-Dm||||m|| + ||-D 2 m|| ||n|| + ||.Dm||||.Dn|| + ||m|| ||-D 2 n|| 
+ IMIII-EH + ||-Dn||||n||} 
< C(l + t)-i(A!(t)) 2 + C(l + _1 (Ai(*)) 2 . (3.63) 

In a similar fashion, we are able to estimate the high order derivatives for density 
as follows. Indeed, in terms of 03.52p . the Holder's inequality and Nirenberg's 
inequality 

< C || Du || ^ || D 2 it 2 

and the fact 

\\D k E\\ < CWD^nW, k > 1 (3.64) 
due to the Riesz potential representation, we can estimate the term Dn as 

\\Dn(t)\\<\\D(N*U )(t)\\+ [ \\D(m(t-T)*Q(U))(r)\\dT 

Jo 

— t 

< C6{1 + t)"t + J* \\D{m{t - t) * Q{U){r))\\dT + J \\D{m{t - r) * Q{U){r))\\dr 

< C6(l + t)~l + C [\l+t- r)- T H\\Q(U)(r)\\ L1 + \\DQ(U)(r)\\)dr 

Jo 

+ C J \\m(t - r) * DQ(U)(r)\\dr 

< CS(1 + t)"t + C(A 1 (t)) 2 [ 2 (1 + t - r)~i (1 + r)- 1 ^ 

+ c Ai + t-r)- 1 !^^)^)!!^ 

"'I 

t 

< CS(1 + + C(Ai(t)) 2 / 2 (1 + t - r)-i (1 + r)" 1 ^ 

Jo 

+ C^t)) 2 J t (l + t- r)~\l + r)-idr 

< C5(l + t)~l + C(l + *)~*~ e (Ai(i)) 2 , (3.65) 

where and below e > is a small but fixed constant. As for D 2 n, it is easy to 
verify that the nonlinear terms in right hand side dominating the time-decay rate 
consists of D(D 2 (mn)), D 2 (Dm-m) and D 2 (D 2 (mn)), which can be estimated due 
to the following facts 

\\D 3 n ■ m(t)\\ < \\D 3 n(t)\\\\m(t)\\ LOO < C5(l + t)~l A^t), 
||£> 3 m -m(t) || < ||D 3 m(t)||||m(t)|| L co < C8{1 + fp*Ai(f), 
||D 4 n-m(t)|| < ||£> 4 n(*)llll™(*)IU» < C5(l + *)"*Ai(*). 
Thus, we can obtain after a straightforward computation that 

\\D 2 n(t)\\ < \\D 2 (N*U )(t)\\+ f \\D 2 {m{t-r)*Q{U)){r)\\dr 
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< C6(l + t)-Z + [ 2 \\D 2 {m{t - r) * Q{U){r))\\dT 

Jo 

+ j\\D\m{t - r) * Q(U)(r))\\dr 

t 

< cS(l + t)-i+C [\l + t -rr'i(\\Q(U)(r)\\ Ll + \\D 2 Q{U){r)\\)dr 
Jo 

\(l+t-rrH\\Q(U(r))\\ + \\D'Q(U(r))\\)dr. 

< C6(l + t)-* + C(A 1 (t)) 2 / 2 (l+t-r)-I(l + r)-frfr 



+ C5A 1 (t) j (l + t-r)-2(l + t)-4d r 

< C6(l + t)~i + C(A!(t)) 2 (l + t)~i + C5A 1 (t)(l + t)"*, (3.66) 



_ 3 



Next, in terms of Proposition 13.11 (13. 46ft . (I3.49[) . (I3.53p . the a-priori assump- 
tion (13.591) . the Holder's and Nirenberg's inequalities, we can prove the time decay 
rate for m and its derivatives as follows. 

\\(m-M*U )(t)\\ < [ t \\M(t-r)*Q(U)(r)\\dT 

Jo 

<C f\l+t-rrH\\Q(U)(r)\\ + \\Q(U)(r)\\ Ll )dr 
Jo 

<C [ (1 + 1 — r)~3(A 1 (t)) 2 (l + r)~ l dr 
Jo 

^Cil+ty^iA^t)) 2 , (3.67) 
where and below e > is a small but fixed constant, and 

||2>m(t)|| < \\D{M * U )(t)\\ + / \\D(M(t - r) * Q{U){r))\\dr 

Jo 

< CS(1 + t)~l + C [ 2 (1 + t - t)-^(\\Q(U)(t)\\li + \\DQ(U)(r)\\)(dr 

Jo 

+ C ^(l + t-r)-«-'(||g(tO(r)|| L i + ||I)g(tO(r)||)dr 

< CS{\ + tyi + C(l + t)-'^ £ (A 1 (t)) 2 , (3.68) 
where we have used ( I3.63p . As for D 2 m, we have 

\\D 2 m{t)\\ < \\D 2 (M*U )(t)\\ + f \\D 2 {M{t-T)*Q{U){T))\\dT 

Jo 

t 

< CS(1 + t)~i + C [° (1 + t - r)-J(||Q(l7)(T)|Ui + p 2 Q(C/)(r)||)dr 

^(IIW)(r)IUi + P 2 g(t/)(r)||)rfr 



+ C j t _ 0- + t-r) 



3_ I 

"4 2 1 
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< C8{\ + t)~i + C(l + t)-f- e (A 1 (t)) 2 + Ctf + 1 - r)~t-3(l + r )-tdr 

< CtfAi(*)(l + t)~t + C(l + i)-t- e (Ai(t)) 2 + C5Ai(t)(l + t)"f . 

(3.69) 

What left is to obtain the time-decay rate for E in terms of (13.641) and (13.471) . 
Indeed, it is easy to get 

\\(E — L* U )(t)\\ < f ||£(t - r) * Q([/)(r)|Mr 
Jo 

< C f\l + t - r)-i(||0(Z7)(r)|| + \\Q{U){r)\\ L ,)dr 
Jo 

<C I (l + t-r)-HA 1 (t)) 2 {l + r)- 1 dr 
Jo 

<C(l + t)-3- £ (A!(t)) 2 . (3.70) 

Step 2: The higher order energy estimates. To enclose the a-priori estimates 
and prove the claim (I3.6ip . we need to derive the time decay rate of (n,m) with 
respect to higher order derivatives as in [2H EH1 [22]. To this end, we look on 
the compressible Navier-Stokes-Poisson system (I3.27l) - (l3.29p as the compressible 
Navier- Stokes with nonlinear inhomogeneous terms related to the electric field E. 
Indeed, by (ll.ip -( Ti~4l) . we derive the system for (p,u) — (1 + n, — ) 

d t n + V-u = f 1 , (3.71) 
d t u + Vn + E - fnAu - /i 2 V(V- u) = -f 2 - / 3 , (3.72) 
E = V(-A) _1 n, lim £^0, (3.73) 

\x\— >oo 

n(x, 0) = n (x) =: p ( a; ) — 1? M ( a; ; 0) = «o(x) = m o( x )/Po( x )i x G M 3 , (3.74) 
where 

/i = u, d x n, d x u) =: -nV- u-Vn-u, 
h = f2{n, m, d x n, d x m, d^n, dim) 

=: u • V u + 1 - y \ > )Vn + /ix — — )Au + /i 2 — — V V- u ), 
1 + n 1 + n 1 + n 

f 3 = -raV$ = -nE, 

and we have chosen p = l,p'(l) — 1, p — pi, (p + v) = p 2 and d — 1, for simplicity. 
Notice that, in view of (13.7ip and (13.721) . there are only two terms different from 
the classical Navier-Stokes equations, that is, E and nE. 

Taking J (13.711) x ndx + f (13.721) • udx and integrating the resulted equation by 
parts, and making use of the facts E = V$ and 

V$ • udx = - <£V • udx = / $n t dx + / <E>/iote 
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$V ■ E t dx + J <$>fxdx = J E ■ E t dx + J Qfrdx 

-lit I \ E \ 2dx -\\^\\fA^lj t J l^l 2 ^-cil^HIIMI L § 

(3.76) 

due to the Holder's inequality and Sobolev inequality, we can obtain after a straight- 
forward computation that 

1 d 
2di 



(|n| 2 + \u\ 2 + \E\ 2 )(t)dx + C J \Du(t)\ 2 dx 

< J fm-ih + fo-udx + CWEWWhW^. (3.77) 
Since it is easy to verify 

\\h\\ L i<C\\{nDu + uDn)\\ Li < C\\n\\ L ,\\Du\\ L 2 + C\\Dn\\ L 4u\\ L s 

< C(\\n\\ L 2 + \\n\\ L 6)\\Du\\ L 2 + C||Dn|| L2 (||n|| L 2 + ||u|| L6 ) 
<C(\\n\\ + \\Dn\\)\\Du\\ + C\\Dn\\\\u\\, (3.78) 

we can control the last term on the right hand side of (13.771) as 



I 2 

< C\\E\\ || (n, Dn, Du) || 2 + C\\Dn\\ \\ (u, E)\\ 2 . (3.79) 
Substituting (I3.79P into (13.771) we have 

lit S (|n|2 + |m|2 + \ E \ 2 ^ dx + cil^WII 2 ^ 

<C|| J E(t)||||(n, J Dr2, J Dn)(t)|| 2 + C|| J Dn||||(n,E)(t)|| 2 + J f in -(f 2 + f 3 ). U dx 
<CA 1 (t)\\(n,Dn,Du)(t)\\ 2 + C(l + t)^A 1 (t)\\(u,E)(t)\\ 2 

+ J fm~ (f 2 + fs) ■ udx. (3.80) 

Due to the facts 

/i = -nV-«-Vn-«, f 2 ~ 0(l)(nDn + uDu + nD 2 u), f 3 = nE, 
it is easy to estimate the last term on the right hand side of (13.801) as 

fin ~ {h + f3) ■ u dx 



„ Lf < C\\E\\(\\n\\ 2 + \\Dn\\ 2 + \\Du\\ 2 ) + C\\Dn\\(\\u\\ 2 + \\E\\< 



^aeHZJuC*) || 2 H- C7Ai(t) || (ra, i^n, || 2 + a-H^Ct) ||« 

<C(e\\Du(t)\\ 2 + A^^D^Du)^ 2 + ^ + t)'^A i mum 2 ), (3-81) 
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which together with (13.801) and the smallness of constant s > gives rise to 

\t I (^l 2 + l M ' 2 + \ E \ 2 )^) dx + C\\Du{t)fdx 
<CAx(0 1| (n, Dn, Du)(t) f + C(l + t)-*k x {t) \\ (u, E)(t) f 

+ C(l + t)-U 1 (t)||n(t)|| 2 . (3.82) 

In a similar procedure, we are able to deal with the higher order derivatives of 
the solution as follows 

4- [ (\D k n\ 2 + \D k u\ 2 + \D k E\ 2 )(t)dx + C\\D k+1 u(t) fdx 



dt . 

<CA 1 (t)(||n(t)||^ + \\Du{t)f Hk ) 

+ C(l + t)-tA 1 (t)||(n,E)(t)|| 2 + C(l + t)-lA 1 (t)||n(t)|| 2 , (3.83) 

with fc = 1, 2, 3, 4, which implies 

|(IK*)||^ + \\u(t)(t)\\ 2 Hi + p(t)(t)||^) + c||zm*)II^ 

+ C(l+t)^A 1 (t)||(n,E)(t)|| 2 + C(l + t)^A 1 (t)||n(t)|| 2 . (3.84) 

Taking inner product between (13.721) and Vn, integrating the resulted equation 
by parts over M 3 , and making use of J E ■ Vndx — — J V ■ Endx = j \n\ 2 dx, we 
have after a straightforward computation that 



— f u ■ Vndx — J uWn t dx + J \S7n\ 2 + |ra| 2 (ix 



d_ 

dt 



<^ + A l m(n,Dn)(t)\\ 2 + C(l + A 1 (t))\\(Du 7 D 2 u)(t)\\ 2 , (3.85) 

which leads to 

«-V«fa + /|Vn| 2 + N^ 

<C J nVn^x + C , (l + A 1 (t))||( J Du, J D 2 n)(t)|| 2 

<C J uV(V ■ u + nV ■ u + uVn)dx + C(l + A 1 (t))\\(Du, D 2 u)(t)\\ 2 

< CA 1 (t)\\(n,Dn)(t)\\ 2 + C(l + A 1 (t))\\Du(t)\\ 2 . (3.86) 

Similarly, taking J D a (\3J2$) ■ VD a ndx with 1 < \a\ < 3, we can have after a tedious 
and complicated calculation that 



D a u ■ VD a ndx + / I VD a n\ 2 + \D a n\ 2 dx 



d_ 

dt 

^AiWUnWH^ + CCl + AiC*))!!^)!!^, (3.87) 
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which together with (I3.86H and the a-priori smallness assumption of Ai(i) gives rise 
to 

j t J D a u ■ VD a ndx + \\n{t)\\ 2 m < G\\Du(t)f m , < |a| < 3. (3.88) 

Taking the summation (I3.84p + (3 x (I3.88P with (3 > small enough, and with 
A(t) rather small we have 

j t K(t) + C(\\n(t)\\^ + \\Dum 2 H *) 
<C(l + t)-*A 1 (t)\\(u,E)(t)\\ 2 + C(l+tylA 1 (t)\\u(t)\\ 2 
<C(l + t)-5Ai(t)/ir(t), (3.89) 



where 



K(t) = (\\n(t)\\ 2 m + \\u(t)(t)\\ 2 m + \\E(t)(t)\\%< +P I ^ D°« ■ VU-ndx) 

l«l<3 

and 

c(||nW||l, 4 +||^)||l, 4 + ||£;(t)|| 2 H4 ) <^)<c"(ll^)ll^ + lkWII^+ll^(OII^)- 

From (13.891) we have by the Gronwall's inequality that 

K(t) < Ceti^r* A ^\\(n ,u ,E )\\ H . < CS. (3.90) 
This together with (13.891) also leads to 

/ \\(n, Du)(r)\\ 2 Hi dT < C6, (3.91) 
Jo 

and finally 

\\(n,m,E){t)\\ H *<C6. (3.92) 

Step 3: Closure of the estimates A3. 61\) . The combination of (13.621) . (13.651) - 
dTTOD and (JS12D leads to 

Ai(t) < C5 + C5A 1 {t) + C{A 1 {t)) 2 , (£[0,3], (3-93) 

which together with the smallness of 5 > leads to the first estimate in (13.611) . 
The proof is completed. □ 

The proof of Theorem 0771 and Theorem \1.2\ The global existence of smooth 
solution of the IVP problem for the compressible Navier-Stokes-Poisson system 
(jl.ip — (jl.4p follows from the short time existence theory, the uniformly a-priori es- 
timates, and the continuity argument. The time-decay rate in Theorem 11.11 follows 
from the Lemma 13.31 The optimal time-decay rate in Theorem 11.21 follows from 
the combination of the Proposition 13.21 and the uniform estimates (13.611) . (13.621) . 
(15m and (137701) . 
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4 L p - time decay rate 
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4.1 L p decay rate for linear semigroup 

In this section, we investigate the L p - time decay rate for the solution of linearized 
NSP system, with p G [2, oo]. To this end, we need to analyze the Green's function 
G = e tB that formed by the linearized NSP system. 

We have the LP time decay rate for the linear semigroup as follows. 

Lemma 4.1 (L p decay rate) Let (n, m) — G * Uq and E = L * Uq with N, M, L 
defined by (13. 38ft and ( 13.431) . then we have 



\D«n{t)\\ LP = \\D«(N*U )(t)\\ LP < C(l + 1)-^ 1 -*^ (\\U \\ L i + \\D«U Q \\ 



\\D«m(t)\\ LP = \\D«(M*U )(t)\\ LP < C(l+ty^ L P )+ ^^(\\U \\ L i + \\D«U \\ LP ), 
\\D a x E{t)\\ LP = \\D:(L*U )(t)\\ LP <C(l+t)--i {1 - l P )+ ^(\\U \\ L i + \\D:u \\ Lp ), 
where p G [2, oo], and the constant C > independent of time. 

To prove Lemma 14.11 we shall make use of following two lemmas [H] . 

Lemma 4.2 (L p multiplier, [8]) Let n > 1 and assume that /(£) G L°° fl 
C n+1 (R n /{0}), with 

\D?f(0\ <C \ , (4.1) 

where ax, o~i > and n > 2 — 2o"2. TTien /(£) zs continuous at both and oo, and 

f = m 1 + m 2 5, 

where ni\ G L 1 (R n ) satisfies \\mi\\ < C(C'), ni2 is the constant 

m 2 = (2tt) - £ lim /(£), 

|£|-+00 

and 5 is the Dirac distribution. In particular, /(£) is a strong L p multiplier, 1 < 
p < oo, in sense that, for any g G LP , 

< C||a|| L p, l<p<oo, 

where C > depends on \m 2 \ < H/Hz, 00 an d iae constant C' . 

Lemma 4.3 ([£]) Lei = K(£,t)f(£) with K(£,t) = e~ c ^ 2t and /(£) G 

L°° n C n+1 (R n ) satisfying 

\Dff(Z)\ <C'|er l/31 , |/3|<n + l. (4.2) 

Then, for any t > D x g(-, t) G L p (M n ) /or a// a and p G [1, oo], and it holds 

ll^("^)IU-(M-) < C(|«|)tt(|"D-f . (4.3) 
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The proof of Lemma 4J_ '■ To make use of the above lemmas to prove Lemma 
14. 1[ we decompose the Fourier transform N, M, L into low frequency term and high 
frequency term below. 

Define 

where (-)i = x(0(')' (O2 = (1 — x(0)(')' an d x(0 is t ne smooth cut off function 
that 

f 1 iei < r, 

X = (4.5) 

\o iei>i2+i. 

Thus, in terms of (13.381) and (14.4H . we have the following decomposition of (n, m) = 
G * Uq in Fourier modes 

n=N-U = N 1 -U + N 2 -U = (M + 9ti) ■ C/ + (A/" 2 + <tt 2 ) • C/ , (4.6) 
m =M X .U + M 2 -U = (Mi + mil) • U + (Af 2 + Sfe) • ^o- (4.7) 

Let we first analyze above higher frequency terms denoted by (-) 2 . Recall that 

A = -/x|£| 2 , (4.8) 

A+ = - (a* + Hl£| 2 + |V4(c 2 |ei 2 + A^)-(2/, + z.)2|e|4 

= -(2^+^)iei 2 + ^; +o(ier 1 ), iei»i, (4.9) 

A_ = - + Ivm 2 - |z v /4( C 2|e| 2 + A- 2 )-(2 Ai + z,)2|e| 4 

= - 2 ^ + o(\^n, iei»i- (4.10) 

We shall prove that the higher frequency terms are L p Fourier multipliers with an 
exponential time decay coefficient Ce~ C5t . For simplicity, we only show that M 2 is 
an LP Fourier multiplier at higher frequency as follows. It holds 

A + e A - t -A_e A + t _ A + e A - t -A-e A - f + (A_e A - i -A_e A + t ) _ \_t , A_e A - t _ A,e A +* 
A+-A_ ~~ A + -A_ ~ e ' A+-A_ A + -A_ ' 

and 



mo = - (a* + f )iei 2 - ^vW+A^vI^m 



- 0* + I )l^l 2 } + i V / (^ + ^) 2 KI 4 -4(c 2 |e| 2 + A- 2 ) 



,2 



*- 2 +(^) 2 

2 ^ | v /( 2 ^) 2 l«l 4 - 4 ( c2 l«l 2 + A " 2 )+(A t +f)l?l 2 -2^7' 

where 



V / f2p + W r W+A 1 ) ~ Odei 2 ), |£| » 1. (4.11) 



Li, Matsumura, & Zhang 

By a direct computation, it is easy to verify 
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|^ a A_| <C\£\ 2+a , |£|»1, (4.12) 

which together with (I4.9p - fl4.10p gives rise to 

f |f | < R, 

l^[(l-x(-))e A -*]| <C{ (4.13) 

' J \e-^|- |a| |£|- 2 \t\>R, 

where and below R > is a given constant, and 

\ At { Ifl < -R, 

l^f 1-X - rV <C< (4.14) 

Thus, from Lemma 14.21 it follows that the inverse Fourier transform of the term 
(1 — x('))( eA_ * + a+-a ) * s an L p multiplier with the coefficient Ce~ C6t . The other 
part of 7V2 at higher frequency can be written as 

(i - = e-<*+* )|€|at ■ ( ^ggg^ a - x(-m 

whose first term in right hand side can be looked on as the function K (£, t) of 
Lemma 14.31 and the rest term satisfies the condition f)4.2l) due to the facts (A + — 
(A*+ f )le| 2 )^ ~ -c 6 |e| 2 ^ for |e| > l and D^e^) < C|£|- fc ((|£| 2 t) + (|£| 2 *) 2 + • • • + 
{\i\H) k )e^\ 2t < C\t\~ k . Thus, the inverse Fourier transform of above term is also 
an LP multiplier with the coefficient Ce~ C(iR2t . These facts imply that A/2 is an LP 
multiplier with the coefficient Ce~ Crt . 

Applying the similar analysis to the terms 9t 2 , M.2-, ^2, £2, and £ 2 , we can 
show that their inverse Fourier transform are all LP multipliers with the constant 
Ce~ C8t . Thus taking C5 = min{c7, eg} and then 

\\D%(N 2 * f),D%(M 2 * /),D»(L 2 * f)\\ LP < Ce- C5t \\D a f\\ LP , (4.15) 

for all \a\ > 0, and p G [2, 00]. 

To prove Lemma 14.11 we also need to deal with the corresponding lower fre- 
quency terms denoted by (-)i. Indeed, we can have 

KMWIU P <c , (i + *)-i( 1 -^, WD^m^KCii + t)-^ 1 '^-^, 

(4.16) 

WD^m^t) \\ LP < c(i + t)-i (1 -l )+ ^-^, < c(i + t)~ f(1 ^^^, 

(4.17) 



(4.18) 



for t > and p G [2, 00] 
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For simplicity, we only prove (14.16P and (14.17p . since the analysis of ( I4.18P 
can be carried out in a similar way. By fl2.17p - fl2.19p . we have the asymptotical 
approximation 

A+eA : t :"- eA+f ~ o(i) e -o*+i")iei a *, A+e * +t :*- eA ~' ~ o(i) e -(M+^)i«i 2 *, i^i « 1, 



A+-A- ^V-V"- i A+-A- 



A + -A_ 



which imply for |£| <C 1 that 

" |M| ~ 0(e~ c ^ 2t ), |oTi| ~ C(e- C9| « |2t |e|), 

1 (4 19) 

|SDTi| ~ C(e^ C9| « |2 '^), IMil ~ C(l)e- C9|?|2f . 

Thus, by Hausdroff- Young's inequality, we can have 

\\D*Afi(t)\\ L r> < C{ [ ||^| |Q| e- C9l « |2 *|^}^ < Cr^-i)-^, (4.20) 

■/|£l<» 



IP^iWIIlp < C{ I ||^| |Q|+1 e- C9l?|2 '|^}^ < Cr*^-*-^, (4.21) 
At\<v 

with - + - = 1 and p G [2, oo], and 

WD^^Wlp <C{ [ W^^e^yd^ <C(1 + *)~ |(1 ~* )+ ^ % (4.22) 
J\i;\<v 

WDZMiWWu, <C{[ ||£|l < V C9 ^ 2t | 9 d£}i < C(l +t)-K 1 -i)-^. (4.23) 

These give rise to lOBjl and lQ7jl . 

Combining the estimates f l4.16l) - fl4.18l) and (14. 15ft . we finally have for t > that 



ra^*/)(<)iii>=ira(M+jv 2 )*/)ii 



<C(l + t)-^-^\\f\\ Ll +Ce- Cet \\D:f\\ LP 

<C(l + f)" i(1 "* )_Jsl (ll/llLi + II^/Hl,), P G [2, oo], 

where JVj = jVj + i = 1,2, and similarly 

\\D2(M * /)(t)||x, < C(l + t)-^-^-^(\\f\\ Ll + ||J3-/||lp), p G [2, oo], 
\\D:(L * f)(t)\\ LP < C(l + t)"i ^->^(\\f\\ L1 + p G [2, oo]. 

The proof of Lemma 14.11 is completed. □ 
In addition to Lemma 14.11 we can also have following estimates, which can be 
shown in a similar procedure as in the proof of Lemma 14.11 and are applicable to 
the L p time decay rate of solution of the original IVP problem for the nonlinear 
Navier-Stokes-Poisson system. Indeed, we have 
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Lemma 4.4 It holds for the Green's functions 9T, M., and £ defined by (13.390 - 
(13T4TD totrf 



pwt */)(*)! 


LP 


< C(l- 




3 /1 


1 \ 

?)- 


i 

2 


I/* 1 


+ P. q /IIlp), 


(4.24) 


ITO* */)(*)! 


LP 


< C(l- 




3/1 


1 \ 


1 

2 






(4.25) 


££CM */)(*)! 


LP 


< C(l- 




3 /1 
-f(l- 


1 S 


<:v 
~2~ 




I^/IIlp), 


(4.26) 


*/)(*) I 


LP 


< C(l- 




3/1 
2 2 


-§)- 


[a] 
2 


(II/IIl 2 + 




(4.27) 


*/)(*)! 


LP 


< C(l- 




-f(l- 


f>- 


a 

2 


(\\f\W + 


I^/IIlp), 


(4.28) 


IIW */)(*)! 


LP 


< c(i- 




3/1 
2 V 2 


-i>- 


M 

2 


(II/IIL2 + 


K/IIlp), 


(4.29) 



with p G [2, oo], \a\ = k > 0, and the constant C > independent of time, so long 
as f e L 2 (R n ) n W k ' p (R n ) or f e L 1 (R n ) D jy fe ' p (M n ). 

Proof: The estimates (fl~2D . (fl~260 and (fl~28l) follow directly from those of 
(I4.16I) - (I4.18I) and the related estimates for higher frequency terms. We need only 
to prove (l4~25l) . (|4~27j) and (l4~29l) . Indeed, it holds for p G (2, oo] that 

\mm*M)\\i* <\\D2Wi * f)h* + * /)Hlp < * f)\\ L v + ce-^\\D:f\\ LP 

and 

\\i%{ni*M)\\z*<c\\(iZ) a vtif\\i* 

<c[{ [ (iei H+1 e^ |€| Y a ^de) 3? -(iii/l , ii,f)]' 

J\e\<v 

<C||/|U2(l + t)-i(H)-i-¥ 
with i + i = 1. As for p = 2, a similar direct computation gives rise to 

* /)(*)IU- < CMWWJWv < C(l + 0"'-^ • II/IIl- (4.30) 

Thus, we obtain the estimate (14.250 . The proof of (14.270 and (14.291) can be com- 
pleted in a similar fashion and the details are omitted. □ 

4.2 LP decay rate for nonlinear system 

In this subsection, we show the IP (p G [2, oo]) time decay estimates for the original 
nonlinear problem. We can verify that the expression of the solution of the IVP 
problem ([3T270 - (!33D]) is 

n = N*U + [ m{t - t) * Q{U){r)dr, (4.31) 
Jo 

m = M*U + [ M{t - t) * Q{U){r)dr, (4.32) 
Jo 
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E = L*U + [ £{t-T)*Q{U)(r)dT. (4.33) 



o 



We have Theorem 11.31 concerned with the LP time decay of strong solutions 
below. 

Lemma 4.5 Under the assumptions of Theorem } l.'J\ the global solutions (n, m, E) 
with E = V$ of the IVP problem (l3~271) -( l3~30l satisfies 

\\n(t)\\ LP < C(l + t)- |(1 "^(||f/ ||Li + \\Uo\\lp), (4.34) 
|| (m, E)(t) Wl, < C(l + t)" f(1 -^ )+ ^(||t/o||Li + \\U \\lp), (4.35) 



with p G [2, oo], Uq = (no, mo), and the constant C > independent of time. 



Proof: By Lemma 13.31 we have for I > 5 that 

\\D k n\\ L 2 < C(l + t)-$-U, k = 0,1,2, \\D 3 n\\ L 2 < C(l + t)~ 1 5, 

\\D k {m,E)\\ L 2 < C(l +t)-*~U, k = 0,1,2, ||L> 3 n|| L 2 < C(l + t)~ 1 5, 

where 5 =: ||(po — P, m o)\\H l (w s )nL 1 (R :i ) > as defined in Theorem 11.31 

Noticing that the nonlinear term Q(U) — V • H contains mainly those like 
0{mD 2 n, Dm ■ Dn, nD 2 m, mDm, nDn, nE}, and using Sobolev embedding theo- 
rem, we can obtain 

\\Q(U)\\ L i ^Cd^l+trWlQWUr, <C5\l+t)- 7 *, pe [2,oo]. (4.36) 
Thus, by (14.311) . (14. 36p . Lemma [4.11 and Lemma [4.41 we are able to show 

\\n(t)\\ LP <\\N*Uo(t)\\ LP + [ \\m(t-r)*Q(U)(r)\\ LP dr 

Jo 

<C(l + ty^ l P \\\Uo\\ L i + \\Uo\\Lv) + C [ \\Af(t-r)*Q(U)(r)\\ LP dT 

Jo 

^(l+tr^aitfoik+ii^oM 

+ c !\i + 1 - Tr^-H\\Q(U)\W + \\Q(U)\\ LP )dr 







3/ 1 _ 1\_ 1 



+ Cj t (l + t-r)--^-p-H\\Q(U)\\ L2 + \\Q(U)\\ LP )dr 

<c(i + tr* (1 -*)(||t/- |Ui + h^oIUp) + c(i + ty^-^s 2 

^Cil + ty^'p^, (4.37) 



and 



/ \\M(t-r)*Q(U)(r)\\ LP dT 
o 
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<C(1 + t)~^ (1 ^^(||f/ ||Li + \\UoWlp) + C [ \\M(t - r) * Q(U)(r)\\ LP dr 

Jo 

<c(i + ty^^(\\u \\ L i + \\u \\ LP ) 

-f(i-l) 



+ C / (l + t- r )-^(||g([/)|| Ll + ||Q(t/)|| iP )dr 
Jo 

+ c J t (i + 1 - t)-s^-?5(ho(z7)IIl» + ng(toilLp)*- 

<C(1 + t)" |(1 "^ )+ ^(||t/ ||Li + II^oIIlp) + C / 2 (1 + t - r)-t (1 -^(l + rY l dT 

Jo 

ft 

+ C (l+t-ry^-p\l + T)-idT 

<c(i + + c(i + ty^ {1 ~p )+ ^- £ 5 2 

<C(l+ty'^p )+ h, (4.38) 

where p £ [2, oo], C > is a constant independent of time, and e > is a small but 
fixed constant. 

The estimates of can be obtained in terms of ||m|| and Riesz potential 

as 

\\E(t)\\ LP < C(l + ty^ 1 -^ + -y\\U \\ L i + \\U \\ LP ), pe [2,oo]. (4.39) 
The proof of the Lemma 14.51 is completed. □ 
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